Correlation Differences in Heartbeat Fluctuations During Rest and Exercise 
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We study the heartbeat activity of healthy individuals at rest and during exercise. We focus on 
correlation properties of the intervals formed by successive peaks in the pulse wave and find signif- 
icant scaling difi'erences between rest and exercise. For exercise the interval series is anticorrelated 
at short time scales and correlated at intermediate time scales, while for rest we observe the oppo- 
site crossover pattern — from strong correlations in the short-time regime to weaker correlations 
at larger scales. We suggest a physiologically motivated stochastic scenario to explain the scaling 
differences between rest and exercise and the observed crossover patterns. 

PACS numbers: 87.19.Hh, 05.45. Tp, 89. 75. Da 

One of the important questions in the analysis of com- 
plex physiological time series is how such series reflect the 
dynamical properties associated with the underlying con- 
trol mechanism 0,^. Recently, it was found, e.g., that 
the fluctuations of the heart interbeat intervals reveal 
long-range power-law correlations ||^ and hidden scale- 
invariant structure Q which may be useful for diagnosis 
and prognosis Here we study the correlation (scal- 
ing) properties of heartbeat dynamics as reflected by the 
pulse wave measured from the finger . 

Previous studies of long interbeat interval series have 
focused primarily on 24h Q and 6h ||^ records, which 
include periods of rest as well as periods of a more in- 
tensive physical activity. However, heartbeat dynamics 
can change dramatically with physical activity. Thus im- 
portant differences in cardiac regulation associated with 
rest and exercise may not be clearly seen when analyz- 
ing records which mix together rest and exercise regimes. 
Here we consider rest and exercise activities separately. 
We focus on the correlations in the interpulse interval 
(/Pi) series derived from the pulse wave signal during rest 
and exercise (Fig. pi). By studying the changes in the 
correlation properties we wish to achieve a better under- 
standing of the physiological mechanism that regulates 
heartbeat dynamics at rest and during physical exercise. 

We analyze 21 records from healthy subjects. Each 
record includes 4 different stages of physical activity de- 
noted as restl, exercisel, rest2 and exercise2 (Fig. ^). 
At the first stage (restl) we measure the /P/ under nor- 
mal rest conditions. At the next stage (exercisel) sub- 
jects are asked to run on a treadmill. After a short re- 
covery, during which subjects sit down to recover their 
heart rate, a new rest-exercise episode (denoted as rest2 
and exercise2) is followed. 
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FIG. 1. (a) A typical example of a pulse wave measured 
as a function of time. As in the case of the electrocardiogram 
signal where interbeat interval fluctuations are studied (see, 
e.g., [^), we analyze the interpulse intervals (IPI) between 
successive peaks in the pulse wave, (b) IPI series obtained 
from the pulse wave signal shown in (a). Each record in- 
cludes two rest and two exercise stages. The duration of each 
stage varies from subject to subject and is between 6-10 min- 
utes, (c) Sign series obtained from the increments AIPI in 
the interpulse intervals during rest and (d) during exercise. 
Note that the sign series of the exercise regime exhibits more 
frequent alternations (stronger anticorrelated behavior) com- 
pared to rest. 

To study the correlation properties of the IPI series 
during rest and exercise stages we use the detrended fluc- 
tuation analysis (DFA) 110] which is a method developed 
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to avoid spurious detection of correlations that are ar- 
tifacts of trends related to nonstationarity. The DFA 
procedure consists of the following steps. We first in- 
tegrate the IPI series to construct the profile Y{k) = 

J2i=i(^P^i ~ {^PI)) where {IPI) is the series average. 
Next, we divide the integrated signal, Y{k), into equal 
windows of size n and find the local trend in each window 
by a least-squares polynomial fit. The order of the poly- 
nomial fit specifies the order of the DFA Then 
we calculate the average of the square distances around 
the local trend. This procedure is repeated to obtain the 
root mean square fluctuation function F{n) for differ- 
ent window sizes n. A power-law relation between F{n) 
and n, F{n) ~ n", indicates the presence of scaling in 
the series. According to random walk theory, the scal- 
ing exponent a is related to the autocorrelation func- 
tion exponent 7 (C(n) ~ when < 7 < 1) and 
to the power spectrum exponent /? [S[f) ^ 1//'^) by 
a = 1 - 7/2 = (/3 -1- l)/2 The value a = 0.5 in- 

dicates that there are no (or finite-range) correlations in 
the data. When a < 0.5 the signal is anti- correlated, 
meaning that large values are most probable to be fol- 
lowed by small values. The case of a > 0.5 indicates the 
existence of persistent behavior in the time series, mean- 
ing that large values are most probable to be followed 
by large values. The higher a is, the stronger are the 
correlations in the signal. 

We also study the correlation properties of the sign 
series sign{AIPI) derived from the IPI increments 
AlPh = IPh+i - IPh Fig. |l|c and Fig. 0d show 

representative examples of sign series obtained from rest 
and exercise stages respectively. For exercise, the signs 
of IPI increments tend to alternate rapidly, indicating 
a strong anticorrelated behavior. At rest stage, on the 
other hand, the signs alternate every several points, and 
thus this dynamics may be characterized by a more cor- 
related behavior. 

Due to the fact that during exercise the IPI series ex- 
hibits strong short-range anticorrelations, we first inte- 
grate the IPI series for all rest and exercise episodes (in 
addition to the integration built in the DFA method), 
to avoid inaccurate estimation of the scaling exponents 
for exercise segments. Because of the apparent linear 
decrease of the IPI during the exercise stage (Fig. |l|b), 
the extra integration introduces a parabolic trend. To 
eliminate the effect of this parabolic trend in the exercise 
stage, we perform 3'^'^ order DFA on the integrated 
IPI series. The integration procedure is not necessary 
for evaluating a for the rest episodes, since they exhibit 
correlated behavior |p^ . 

In Fig. Ha we present the fluctuation function, F{n), 
of the integrated IPI series for rest and exercise seg- 
ments of a typical subject. For all 21 individuals we 
observe a characteristic crossover around n w 20 where 
there is a change in the correlation behavior between 



short and intermediate scales regimes. We denote the 
scaling exponent of the short-range regime as ai (esti- 
mated for scales 8 < n < 14) and the scaling expo- 
nent of the intermediate regime as a2 (estimated for 
30 < n < 300). The type of crossover is different for 
rest and exercise: for the rest ai > a2, while for the 
exercise ai < a2- The fluctuation functions for rest and 
exercise stages construct a "fish" -like curve (Fig. 0). 
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FIG. 2. Fluctuation function, F(n), as a function of time 
scale n (in beat number) for rest (A) and exercise (o) stages 
of a typical healthy subject for (a) the original /P/ series, and 
(b) the sign series sign{AIPI) . For all records we observe a 
crossover between two different regimes of correlations. The 
dashed lines indicate the boundaries of these regimes in which 
short-range scahng exponents ai and intermediate exponents 
Qf2 have been calculated. Note the different crossover patterns 
for rest (qi > Q2) and exercise (ai < C12) stages. 

We apply a similar scaling analysis to sign series de- 
rived from rest and exercise segments of the IPI signal 
(Fig. IJd). The sign series do not have any global trend, 
like the original IPI series have; thus, in this case it is 
enough to use 2"^^ order DFA. For sign series we cal- 
culate the short-range scaling exponent ai in the range 
6 < n < 12 and intermediate exponent a2 in the range 
20 < n < 60 [|0|. 
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-0.14 ± 0.17 
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-0.21 ± 0.06 


0.41 ± 0.1 
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1.21 ± 0.25 


0.91 ± 0.12 


0.23 ± 0.15 


0.22 ± 0.07 



TABLE I. Comparison of scaling exponents Qi and a2 
between rest stages, exercise stages and whole records that 
include rest and exercise episodes altogether. For each stage 
the average scaling exponent ± standard deviation are shown. 
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We obtain a complete separation between rest and 
exercise for the original and sign series in the short 
range (Fig. ||) and an almost complete separation for the 
intermediate-range scaling exponents (Table |). The p- 
values for the original and sign series (obtained by the 
paired samples Student's t-test are less than 10^^" 
for the short-range regime and less than 10"'' for the in- 
termediate regime. We find that our results are robust 
and do not change significantly with repetitive rest and 
exercise stages (Table I). 
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FIG. 3. Short-range scaling exponents a\ from 21 healthy 
individuals (a) for the original IPI series and (b) for the sign 
series at rest (•) and during exercise (o). At the right hand 
side we show the average ai ± standard deviation. In accor- 
dance with Fig. ^ the short-range exponents during exercise 
are significantly smaller than during rest. Note the complete 
separation of the two stages, emphasized by the dashed lines. 

To illustrate the importance of considering separately 
rest and exercise episodes, we perform our analysis also 
on the entire IPI records which include rest and exercise 
episodes altogether (Fig. |l|b). Wc find indeed that the 
scaling of the whole record reflects neither the correlation 
properties of rest, nor of exercise (see Table |). 

The significant differences between the values of ai 
(Fig. H) and the different crossover patterns for rest and 
exercise stages (Fig. ^) may offer some insight on the un- 
derlying physiological mechanism controlling the heart- 
beat dynamics. The heart rhythm is regulated mainly 
by the parasympathetic (PS) and the sympathetic (SM) 
branches of the autonomic nervous systems pl[ |. PS im- 
pulses slow the heart rate while SM impulses accelerate it. 
The interaction between these two branches is reflected 
by the time organization of the IP/ series (Fig. |l|b). 

The principle of homeostasis (dynamic equilibrium) as- 
serts that physiological systems seek to maintain a con- 
stant output in spite of continuous perturbations [p2| . 
However, healthy systems even at rest display highly ir- 
regular dynamics (see Fig. lb) 0J|,|2^. In a recent work 
Ivanov et al. proposed a general approach based on 



the concept of stochastic feedback to account for the 
complex fractal variability in biological rhythms. In this 
framework the time evolution of a physiologic system, 
e.g. the heartbeat dynamics, can be represented by a 
random walk biased toward some preferred "attracting" 
levels. Both the SM and PS systems controlling the heart 
rhythm generate attracting levels which bias the walker 
(modelling the interbeat interval series) in opposite di- 
rections leading to complex heartrate fluctuations. Al- 
though these attracting levels change in time, according 
to the response of the intrinsic physiological mechanism, 
they can vary in a limited range only, thus keeping the 
walker away from extreme values. 
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FIG. 4. Schematic illustration of a random walk with (a) 
two different levels of attraction and (b) a single attraction 
level, and their sign decomposition (c) and (d). (a) The com- 
bined effect of the upper level (representing the "preferred" 
level of parasympathetic (PS) system) and the lower level 
(sympathetic (SM) system) do not restrict the walker's fluc- 
tuations for the short-time regime. But for the larger times 
the bounded walk results in a crossover to a less correlated 
behavior, similar to the crossover obtained for the rest stage 
(Fig. (b) On the other hand, the walker attracted by a 
single SM level produces short-range anticorrelations with a 
crossover to a more correlated behavior in the intermediate 
regime due to the alterations of the SM level (mimicking ex- 
ercise) . The sign series for the two attractive levels presented 
in (c) is more likely to cluster than the series for a single level 
scenario presented in (d) (Cornpare with sign dynamics for 
rest and exercise shown in Fig. hlc and Fig. hid). 

Based on this general approach we suggest a schematic 
scenario that explains the different crossover patterns of 
the /P/ fluctuations for rest and exercise (Fig. |^). At rest 
both the SM and PS systems are active, and each of them 
attracts the walker toward its own level (Fig. 4a). Since 
the response time of the PS system is shorter than that 
of the SM system ||2^ , ^ , we assume that the preferred 
attracting level of the PS system alters more rapidly 
that the one related to the SM system. This scenario 
can account for the crossover observed for the rest stage 
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(Fig. 2). When the walker is between the two attracting 
levels, each level imposes a bias in an opposite direction. 
Thus the walker is able to move in both directions until 
he crosses any of the two levels after which he is pulled 
back. This model scheme reproduces the crossover in the 
scaling behavior (Fig. 2) from a larger value of the cor- 
relation exponent at short scales, where the fluctuations 
of the walker are not bounded, to a lower value of the 
exponent at large time scales, where the dynamics of the 
walker is bounded by the SM and PS attracting levels. 

During exercise the SM system dominates p6| and 
the dynamics can be described effectively by a single 
attracting level (Fig. 4b). In this case the walker fluc- 
tuates around this level producing an anticorrelated be- 
havior at short time scales. However, since the attract- 
ing level changes with time and since the walker follows 
these changes, the fluctuations in the random walk in- 
crease at intermediate time scales, causing a crossover to 
a more correlated behaviour. This scheme accounts for 
the observed crossover pattern in the scaling of the IPI 
fluctuations from an anticorrelated behavior with small 
value of the corretalion exponent at short time scales to 
a correlated behavior characterized by a larger value of 
the exponent at large time scales (Fig. 2). Our scenario 
can also explain the remarkable difference in the ampli- 
tude of the fluctuations at rest and during exercise (see 
Fig. lb). When two attracting levels bias the walker (a 
situation in our scenario corresponding to rest) the fluc- 
tuations are larger compared to the exercise stage when 
there is a single dominant attracting level (Fig. 4). 

A "fish" structure similar to Fig. |^ (but with different 
scaling exponents) was observed when comparing healthy 
subjects with congestive heart failure patients Q . These 
results support our scenario, since for heart failure pa- 
tients there are evidences of SM dominance [ p7| , resem- 
bling the state of the autonomic nervous system under 
physical exercise. 

In summary, we study correlations in heartbeat fluc- 
tuations during rest and exercise. We show that the 
significant scaling differences and the different crossover 
patterns between rest and exercise (Fig. |[) can be ex- 
plained based on the "attractive levels" scenario. We, 
therefore, conclude that the interaction between the com- 
peting branches of the autonomic nervous system under- 
lies the correlation properties of heartbeat. 
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